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0\ . Abstract 



We consider Levy flights characterized by the step index / in a quenched 

random force field. By means of a dynamic renormalization group analysis 

^ I we find that the dynamic exponent z for / < 2 locks onto /, independent of 

O ■ dimension and independent of the presence of weak quenched disorder. The 

critical dimension, however, depends on the step index / for / < 2 and is 

given by d^ = 2f — 2. For d < dc the disorder is relevant, corresponding to a 

/\f • non trivial fixed point for the force correlation function. 
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There is a current interest in the dynamics of fluctuating manifolds in quenched 
random environments [0. This fundamental issue in modern condensed matter 
physics is encountered in problems as diverse as vortex motion in high temperature 
superconductors, moving interfaces in porous media, and random field magnets and 
spin glasses. 

The simplest case is that of a random walker in a random environment, correspond- 
ing to a zero dimensional fiuctuating manifold. This problem has been treated 
extensively in the literature 0, |^, ^ and many results are known. 

In the case of ordinary Brownian motion, characterized by a finite mean square step, 
in a pure environment without disorder, the central limit theorem [Q implies that 
the statistics of the walk is given by a Gaussian distribution with a mean square 
deviation proportional to the number of steps or, equivalently, the elapsed time, i.e., 
the mean square displacement 

(r'it)) ex Dt'/\ (1) 

where the dynamic exponent z = 2 for Brownian walk and D is the diffusion coeffi- 
cient. 

There are, however, many interesting processes in nature which are characterized 
by an anomalous diffusion with z ^ 2 due to the statistical properties of the envi- 
ronments [|, ^. Examples are found in chaotic systems 0, turbulence 0, §, fiow 
in fractal geometries 0, and Levy fiights 0, |Tl|], which generally lead to enhanced 



diffusion or superdiffusion with z < 2. We note that the ballistic case corresponds 
to z = 1. The other case of subdiffusion or dispersive behaviour with z > 2 is 
encountered in various constrained systems like doped crystals, glasses or fractals 



g, o, pi ra, IT 



Irrespective of the spatial dimension d, ordinary Brownian motion traces out a 
manifold of fractal dimension dp = 2 |1^. In the presence of a quenched disordered 



force field in d dimensions the Brownian walk is unaffected for d > dp, i.e., for d 
larger than the critical dimension dc = dp the walk is transparent and the dynamic 
exponent z locks onto the value 2 for the pure case. Below the critical dimension 
dc = 2 the long time characteristics of the walk is changed to subdiffusive behaviour 



with z > 2 |jT8|, [T|, 13]- In (i = 1, {r (t)) ex [logt]"^, independent of the strength of 



the quenched disorder [El 



Levy flights constitute an interesting generalization of ordinary Brownian walks. 



Here the step size is drawn from a Levy distribution characterized by the step index 
/ p. The Levy distribution has a long range algebraic tail corresponding to large 
but infrequent steps, so-called rare events. This step distribution has the interesting 
property that the central limit theorem does not hold in its usual form. For / > 2 
the second moment or mean square deviation of the step distribution is finite, the 
central limit theorem holds and the dynamic exponent z for the Levy walk locks 
onto 2, corresponding to ordinary diffusive behaviour; however, for / < 2 the mean 
square step deviation diverges, the rare large step events prevail and determine the 
long time behaviour, and the the dynamic exponent z depends on the microscopic 
step index / according to z = f (/ < 2), indicating anomalous enhanced diffusion. 



that is superdiffusion |11, 22, |2^. The 'built in' super diffusive characteristics of 
Levy flights have been used to model a variety of physical processes such as self 
diffusion in micelle systems |]2^, and transport in heterogeneous rocks ^2 . 



In the present letter we consider Levy flights in the presence of a quenched random 
force field and examine the interplay between the 'built in' superdiffusive behaviour 
of the Levy flights and the pinning effect of the random environment generally 
leading to subdiffusive behaviour. Generalizing the discussion in refs. [p!8| , pO| , p5| we 



find that in the case of enhanced diffusion for / < 2 the dynamic exponent z = f, 
independent of the presence of weak disorder. On the other hand, we can still 
identify a critical dimension dc = 2f — 2, depending on the step index / for / < 2. 
Below dc the quenched disorder becomes relevant as indicated by the emergence of 
a non trivial fixed point in the renormalization group analysis. Here we give a brief 
outline of our arguments while details will be published elsewhere. 

It is convenient to discuss Levy flights in terms of a Langevin equation with 'power 
law' noise [^]. In an arbitrary drift force field F(r), representing the quenched 



disordered environment, the equation takes the form 

^=F(r(t)) + r7(t). (2) 

Here r] is the instantly correlated power law white noise with distribution 

p(r})d'^r] oc r]^^^^dr]. (3) 

We assume an isotropic distribution characterized by the step index /. In order 
to ensure normalizability we have introduced a lower cutoff 77 ~ a of the order of 
a microscipic length a and chosen / > 0. For / > 2 the second moment, (77^) = 
j p{rf)rfd'^ri, is finite and a characteristic step size is given by the root mean square 



deviation J (jf)- For 1 < / < 2 the second moment diverges but the mean step. 



{if), is finite. In the interval < / < 1 the first moment diverges and even a mean 
step size is not defined. 

The power law noise 77, describing the consecutive Levy steps, drives the position r 
of the walker. For the quenched force field F(r) we assume a Gaussian distribution, 

p(F) oc exp (-- / rf'^rrfVF"(r)A"^(r - r')"^F^(r')) (4) 

(F"(r)F'^(r'))F = A"'^(r - r'). (5) 

Here A"^(r — r') is the force correlation function expressing the range and vector 
nature of F(r). 

In the absence of the force field, i.e., F(r) = 0, the distribution for the position 
of the walker is easily inferred. From the definition P(r, t) =< 5{t — r(t)) >, the 
solution of Eq. (2), and averaging according to Eq. (3), we find the scaling form 

1111,113, 



cfik 



P{r, t) = j ^^ exp (zkr - DP | t |) =| t p^ G(r/ | t |^). (6) 

D is a diffusion coefficient setting the time scale. This procedure is equivalent to 
applying the central limit theorem to a sum of random variables with the Levy 
distribution in Eq. (3) [^, |17 . 



The scaling exponent /x depends on the step index /, characterizing the Levy distri- 
bution. For / > 2, i.e., the case of a finite mean square step, yU locks onto the value 
2 and the scaling function G takes the Gaussian form for ordinary Brownian walk, 
G{x) = exp (— x^). This is a consequence of the central limit theorem which here 
leads to universal behaviour. For / < 2 the scaling exponent /i = / and the scaling 
function G can only be given explicitly in terms of known functions for /i = 1, the 
ballistic case, where we find the Cauchy distribution F{x) = {1 + x'^)~^^'^~^^^^'^\ It is, 
however, easy to show that G — > const, for x — >■ and G — > for x — > 00. From the 
distribution in Eq. (6) we deduce the scaling form for the mean square displacement 
of the walker in time t: 

< r^(t) >= f P(r, tyd'^r oc t^ (7) 

and following from Eq. (1) the dynamic exponent z = fi. We also note that the 
fractal dimension of a Levy fiight is dp = fJ^ [ITTf - 



In the presence of the quenched force field given by Eqs. (4) and (5) it is convenient 
to recast the problem given by the Langevin equation (2) in terms of the associated 
Fokker-Planck equation, 

^^^ = -V(F(r)P(r, t)) + DV^P(r, t). (8) 

Here the first term on the right hand side of Eq. (8) is the usual drift term due to 
the motion of the walker in the force field, the second term follows from Eq. (6) 
and the assumption of independent contributions to the probability current. The 
'fractional' gradient operator V^ is the Fourier transform of —k'^ and is a spatially 
non local integral operator reflecting the long range Levy steps; for /i = 2 it reduces 
to the usual Laplace operator describing ordinary diffusion. 

The Fokker-Planck equation in Eq. (8) together with the distribution in Eqs. (4) 
and (5) defines the problem of a random Levy walker in a quenched random force 
field. The 'microscopic' Levy steps occuring on a fast time scale represented by 
the noise term in the Langevin equation (2) has now been absorbed and replaced 
by the anomalous diffusion term in the Fokker-Planck equation (8). The remaining 
randomness due to the quenched force field is assumed static compared to the time 
scale of r(t). 

There are a variety of techniques available in order to treat the random Fokker- 
Planck equation (8). Applying the Martin-Siggia-Rose formalism in functional form 
P7| , P^l , P^ , Pm and using either the replica method [Q or an explicit causal time 
dependence |1^, ^ ^ , one can average over the quenched force field and construct 
an effective field theory. A more direct method, which we shall adhere to in the 
present discussion, amounts to an expansion of the Fokker-Planck equation (8) in 
powers of the force field and an average over products of F(r) according to the 
distribution in Eqs. (4) and (5) 



In order to deduce the scaling properties of the force averaged distribution 
{P{r,t))F and the mean square displacement ((^^(t)))^ "we carry out a renormal- 
ization group analysis of Eq. (8) following the momentum shell integration method 



propounded in refs. [^ 0, g^ and also drawing on the results in ref . |T8[ . 
Defining the Fourier transform 

P(k,cj)= f d'^rdtexpiiujt - ikr)P{r,t)e{t), (9) 



where 6{t) is the step function, we obtain the Fokker-Planck equation 

(_,^ + Dk'^)P{k, uj) = Poik) -ikj 70jF(k - p)P(p, u). (10) 

Here the force field is averaged according to Eqs. (4) and (5), or 

(F"(k)F^(p))^ = A"^(k)(27r)^5(k + p) (11) 

for all pairwise force contractions (Wick's theorem for the Gaussian distribution); 
Po(k) = P(t = 0) is the initial distribution. For simplicity we consider the case 
of isotropic zero range force correlations, i.e., A"^(k) = A^"^^; the general case 
has been discussed in refs.[H, ^, ^. We introduce a microscopic UV cut off and 
assume < /c,p < 1. Iterating Eq. (11), identifying self energy and force correlation 
corrections to second order in A, and performing a momentum shell integration, 
i.e., averaging over the force in the shell e~^ < A;,p < 1, we obtain the 'corrected' 
Fokker-Planck equation 

{-lu + DF + 5De)Pik, u) = Poik) -ikj 7^F(k - p)P(p, u). (12) 

and the force correlation function 

(P"(k)P'3(p))^ = (A + (5A)5°^(27r)'^5(k + p) (13) 

for < A;,p < e~^. For small values of the scale parameter I the corrections 5D and 
5A arising from the momentum shell integrations are proportional to L From the 
diagrammatic contributions to 5D and 5/S. given in refs. |T8|, ^, involving three 
propagators {—iuj+Dk^^)~^ and two force contractions A for 5D and two propagators 
and two force contractions for 5 A, all evaluated in the static limit a; = and on the 
shell fc = 1, we obtain 

A2 A2 

5D = -A—i 6A = -B—i, (14) 

where A and B are geometric factors associated with the area of a d-dimensional unit 
sphere. In order to derive the 'renormalized' Fokker-Planck equation we introduce 
scaled quantities k' = ke^, p' = pe^, cu' = uje°'^^\ P'(k',ct;') = P(k, Ci;)e~"*^^\ and 
F'(k') = F(k)e-('^+^)^+"W such that Q<k',p' <l and find 

^o(k') - ^k' / ^F'(k' - p')P'(p', u') (15) 



and 

(F'"(k')F'^(p'))F = (A + (5A)e-('^+2)^+2"W<5"^(27r)'^(5'^(k' + p'). (16) 

From the renormalized Fokker-Planck equation and force correlation function we 
read off the renormahzation group equations for D and A, 

D' = Z}e^^+"(^) (17) 

A' = (A + 5A)e-('^+2)^+2aW (^g) 

or in differential form setting a{i) = J^ z{i')di' and introducing Eq. (14) 



dDji) 
di 



[z{i) - fi)D{i) (19) 



Before we discuss the renormahzation group equations for the scale dependent diffu- 
sion coefficient D{£) and force correlation A{i) we notice immediately that for Levy 
flights with fi = f < 2 there is no correction to the anomalous term Dk'^ in the 
Fokker Planck equation. The perturbative contribution comes with a leading power 
/c^ which is subdominant compared with the Levy term k^. In the case of ordinary 
Brownian motion for n = 2 there is a correction and we recover the results in refs. 



m. 



Proceeding with the discussion of Eqs. (19) and (20) we fix as usual the dynamic 
exponent z{i) = /i such that the diffusion coefficient stays constant under renormal- 
ization, i.e., D{i) = D. The equation for A(£) then takes the form 



dA{i) ^^ ^ ^^^^^^ ^A(£)2 

15 



{2^i-d-2)A{i)-B-^. (21) 



For d greater that the critical dimension dc = 2^ — 2 Eq. (21) has the trivial 
fixed point A* = and the quenched disorder is irrelevant. For d < dc 'we obtain 
the non trivial fixed point A* = {dc — d)D'^/B and the quenched disorder becomes 
relevant. We also note that unlike the case of Brownian motion the critical dimension 
dc = 2fi — 2 is less than the fractal dimension dp = /i. 



In order to derive the scaling properties of the force averaged distribution 
{P(k,u!))F and the means square displacement {{r'^{t)))F we use the methods dis- 



cussed in refs. [^, ^, |3J]. From the derivation of the renormalization group equa- 



tions we infer the scaling relation 

{P{k,u;,A))F = e'^W(P(ke^a;e"W, A(£)))^. (22) 

In the vicinity of either the trivial fixed point A* = for d > dc = 2fi — 2 or the 
fixed point A* = {dc — d)D'^/B for d < dc ^e have, setting a{i) oc fi£, 

(P(k, uj, A))f = e^^Pike^ cje-^^ A*))^^ (23) 

Choosing ke^ ~ 1 we obtain the scaling form 

{P{k,u;,A))F = k~^L{k/u;'/n, (24) 

where L is a scaling function. From Eq. (24) follows directly 

{{r\t)))F oc t^/^ = f/\ (25) 

Since the scaling analysis in the Levy case is quite similar to the Brownian case 
discussed in refs. ||18|, |^, ^ we shall not repeat this analysis here but simply 
discuss and summarize our results below. 

For Levy flights in random quenched environment we have demonstrated that the 
dynamic exponent z locks onto the scaling index fi, depending on the Levy step 
index /, independent of the presence of weak quenched disorder. The long range 
super diffusive behaviour characteristic of Levy flights enables the walker to escape 
the inhomogeneous pinning environment and the long time behaviour is the same 
as in the pure case. This result is in contrast to the case of a Brownian walk where 
the dynamic exponent z > 2 below d = 2, corresponding to subdiffusive behaviour. 
We have also identified a critical dimension dc = 2^ — 2, depending on the scaling 
exponent /i. Below the critical dimension the weak disorder becomes relevant as 
shown by the emergence of a non trivial fixed point. 

Bouchaud |3l|] has given a heuristic argument yielding the critical dimension dc = fi 



in the Levy case. This result is at variance with the critical dimension dc = 2^ — 2 
given here based on a renormalization group analysis. It would be of interest to 
construe a qualitative heuristic argument for the critical dimension dc = 2fi — 2 
given here and the insensitivity of the dynamic exponent z = fi to the weak quenched 
disorder. 

8 



Below we have in Figure 1 plotted the scaling exponent fi and the dynamic exponent 
z as function of the Levy step index /. In Figure 2 we have plotted the critical 
dimension dc as a function of /i. 



/i,Z {Z = yU,) 



2 - 



Fig.l. Plot of the scaling index /i and the dynamic exponent z as functions of the 
Levy step index /. For / > 2 we have normal Brownian diffusion and we obtain 
corrections to z below d = 2; for < / < 2 we have anomalous Levy superdiffusion. 
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Fig. 2. Plot of the critical dimension dc as a function of the scaling index yU. For 
/i = 2 we have the Brownian case dc = 2; for 1 < fi < 2 dc depends linearly on /x. 
Note that dc = in the ballistic case for /i = 1. 
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